Impurity transport through a strongly interacting bosonic quantum gas 
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Using near-exact numerical simulations we study the propagation of an impurity through a one- 
dimensional Bose lattice gas for varying bosonic interaction strengths and filling factors at zero 
temperature. The impurity is coupled to the Bose gas and confined to a separate tilted lattice. 
The precise nature of the transport of the impurity is specific to the excitation spectrum of the 
Bose gas which allows one to measure properties of the Bose gas non-destructively, in principle, by 
observing the impurity; here we focus on the spatial and momentum distributions of the impurity 
as well as its reduced density matrix. For instance we show it is possible to determine whether 
the Bose gas is commensurately filled as well as the bandwidth and gap in its excitation spectrum. 
Moreover, we show that the impurity acts as a witness to the cross-over of its environment from 
the weakly to the strongly interacting regime, i.e., from a superfiuid to a Mott insulator or Tonks- 
Girardeau lattice gas and the effects on the impurity in both of these strongly-interacting regimes 
are clearly distinguishable. Finally, we find that the spatial coherence of the impurity is related to its 
propagation through the Bose gas, giving an experimentally controllable example of noise-enhanced 
quantum transport. 

PACS numbers: 67.85.-d, 05.60.Gg, 47.60.-i, 72.10.-d 



I. INTRODUCTION 



The transport of an impurity through a Bose gas has 
received much attention due to both its inherent appeal 
as an example of nonequilibrium quantum dynamics and 
its importance in mimicking the transport of an elec- 
tron in a conductor The advantage of a cold atom 
and optical lattice setup is its ability to realize idealized 
Hamiltonians and fine-tune interactions over large pa- 
rameter regimes, often beyond those found in condensed 
matter systems 0. In this vein, we extend the large 
body of theoretical work on the transport of an impurity 

tartly subject to static forcing) through a superfiuid [3|- 
by considering strongly-interacting bosons in one di- 
mension up to the Mott insulator and Tonks- Girardeau 
limits, whereupon the bosons exhibit fermionic charac- 
teristics (Toj . Recently, impurity transport in this con- 
text has received interest due to its successful realization 
by two experiments, one using a species specific dipole 
potential [ll| and another u sing gr avity to provide a 
static force for the impurities jl2l . 1 1 31 ] . In this article we 
look at how an impurity subjected to a one-dimensional 
tilted optical lattice potential moves through a bosonic 
lattice gas in diverse regimes. We simulate the system 
non-perturbatively and near-exactly at zero temperature 
using the matrix product based time-evolution block dec- 
imation (TEBD) algorithm [!l[l5j|. 

This work also relates to progress made in probing 
cold atoms in optical lattices. For example, the exci- 
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tation properties of a Bose gas in an optical lattice can 
be analyzed using Bragg spectroscopy via two pho ton 
scattering [l|| and lattice depth modulation (l7l Il8| - 
the experiment in fl7| having been later successfully sim- 
ulated in [lj| Ho[. Studying the transport of a Bose 
gas after imparting a momentum on it using a magnetic 
field [l8| or tilting the lattice 21 1 has allowed experimen- 
talists to successfully reveal signatures of phase transi- 
tions or determine excitation spectra, respectively. Re- 
cently an impurity has been used in experiments as a 
coherent probe for large quantum systems in order to 
extract impurity-boson collision parameters [22| and in- 
teraction strengths [23| . 

We here combine and extend these two ideas of probing 
a system using impurities and transport; we discuss how 
the motion of an impurity depends on its environment 
and then using our simulations we show what information 
can be revealed about the environment non-destructively 
by observing the impurity. As our example, we consider 
the environment to be a Bose lattice gas and demonstrate 
that by analyzing the transport of the impurity one can 
determine the bandwidth and the gap in the excitation 
spectrum of the Bose gas. In the case of a commensu- 
rately filled Bose gas the impurity propagation provides 
a clear signature of the superfiuid to Mott insulator tran- 
sition in the form of a sharp quench of the expected cur- 
rent. The strong dependence on the filling allows one 
to determine non-destructively whether the bosons have 
been prepared in a commensurately filled state, as desired 
in some applications of quantum information processing. 

The structure of this article is as follows: In Sec. |TT] 
we develop the theoretical framework underlying the dy- 
namics of impurities in a tilted lattice, first decoupled and 



2 




FIG. 1: (Color online) Schematic diagram for the (a) impurity, (b) Bose gas and (c) interaction Hamiltonians. 



then coupled via interactions to an environment. Subse- 
quently, we specialize to the case of an environment con- 
stituting a Bose lattice gas. In Sec. IIIII we simulate the 
propagation of the impurity. Specifically, in Sees. IIII Al 
and IIII Bl we analyze the displacement of the impurity 
in an incommensurately and commensurately filled Bose 
gas, respectively, and show that deviations from the basic 
Esaki-Tsu dependence on tilt clearly reveal information 
about the gas. In Sec. IIII Cl we show how the interaction 
of the impurity with the Bose gas affects the momentum 
distribution of the impurity. Then in Sec. IIII Dl we dis- 
cuss the decoherence of the impurity as an example of 
noise-enhanced transport. We make our conclusions in 
Sec. IIVI and leave a brief review of the TEBD algorithm 
to Appendix lAl 



II. THEORETICAL FRAMEWORK 

A. Impurities in a tilted lattice 

In our system we have one or a few impurities confined 
to a tilted optical lattice which allows them to move along 
one spatial dimension only, depicted in Fig. Ufa). Such 
an optical lattice can be created for the impurities by 
subjecting them to both an optical lattice and a static 
forcing, e.g., gravity fl2| . or by chirping the frequency 
difference between the counter-propagating lasers mak- 
ing up the standin g wa ve of an optical lattice along the 
direction of motion [24j . For sufficiently deep lattices and 
low temperatures solely the lowest Bloch band is occu- 
pied and only tunneling between nearest neighbor sites 
must be considered. The Hamiltonian of the impurities 
is of the form 



H a — —j a 



(1) 



where denotes nearest neighbor sites i and j, a] (a^) 
is the creation (annihilation) operator for an impurity in 
a Wannier state [25[ localized at site i, separated in en- 
ergy and distance from its neighbors by the Bloch energy 
A and the lattice constant a. An additional condition for 
the above Hamiltonian to hold is that the separation of 
the Bloch bands must be much greater than the Bloch 
energy Note that the impurities could be either bosonic 
or fermionic; we assume them to be sufficiently dilute so 



that impurity-impurity interactions and their quantum 
statistics have negligible effects. 

Without a further interaction with an environment it 
is not possible for an impurity to dissipate energy and 
therefore alter its expected position. The eigenstates of 
the Hamiltonian in Eq. ([1} are centered at each of the 
lattice sites i and have energies Ai forming the so-called 
Wannicr-Stark ladder of states. Introducing the width 
A = 2J a /A the creation operator for a Wannier-Stark 
state is defined by d\ = V J}_i(A)aj where J n {h) is a 
Bessel function of the first kind of order n. As the states 
are separated by an energy A a superposition of them 
oscillates at a frequency A/ft. 

This Wannier-Stark picture is consistent with the so- 
lution of the semi-classical equations for an impurity in a 
lattice described by a wave-packet of well defined quasi- 
momentum (therefore spatially it must be spread over 
many lattice sites) [26| . In such a framework the effect of 
the tilt is a constant drift of the quasi-momentum of each 
of the impurities at rate k = A/ ha. Together with the 
semi-classical equation of motion for the group velocity 
of a wave-packet 



IdE 



(2) 



and the single particle energy spectrum in a lattice 
E(k) = 2J a (l — coska) this results in a group velocity 
that is sinusoidal in time. Thus an impurity undergoes 
harmonic motion with frequency A/ h but there is no net 
drift down the lattice. Although originally derived semi- 
classically, these Bloch oscillations are a purely coherent 
phenomenon, can be arrived at quantum-mechanically 
and are rigorously connected to the Wannier-Stark pic- 
In Fig. [5Ja) we show the coherent Bloch oscil- 
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ture 

lation of an impurity initially located in a Wannier state. 

We note that Bloch oscillations in semi-conductors are 
usually obscured by large scattering rates, though they 
have been observed in superlattices where the Bloch fre- 
quencies A/H achievable are much higher due to larger 
lattice spacings (28|. Bloch oscillations are also observ- 
able in analogous optical lattice setups, as we find here, 
due to the much smaller scattering rates. 



B. Impurities coupled to an environment 

If the impurities are coupled to an environment then 
this introduces a mechanism through which an impurity 
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can dissipate energy and fall down the lattice. Esaki 
and Tsu [23 ], using the semi-classical Bloch oscillation 
picture, derived an expression for the average current in 
the relaxation-time approximation. The main assump- 
tion is that an impurity suffers collisions at some rate 1 /r 
that resets its momentum distribution to that of thermal 
equilibrium (we assume zero-temperature here). The re- 
sulting expression for the average velocity of each of the 
impurities is then 



(a) 



(v) 



2aJ a 



-A/h 



h {rA/hf + r 



(3) 



describing a linear (ohmic) dependence on forcing at 
small rA followed by negative differential conductance 
as (rA) -1 when rA ^> 1. This dependence is plotted in 
Fig. Hlb). In this semi-classical framework the negative 
differential current arises if the timescale of the collisions 
is considerably larger than the period of Bloch oscilla- 
tions, in which case the velocity of the impurity is not al- 
ways in the same direction and averages out between col- 
lisions. Negative differential conductance, as with Bloch 
oscillations, is usually obscured in semiconductors but 
due to the much larger values of rA achievable it has 
been observed in optical lattices [30[ and superlattices 

B 

Clearly the relaxation-time approximation is not al- 
ways valid and the actual current will depend on the de- 
tails of the environment and its interaction with the im- 
purity. However, it has been found that at least qualita- 
tively an Esaki-Tsu-like dependence of current on forcing 
is obeyed in many cases: for density-density interactions 
between an impurity and a Bose lattice gas at infinite 
temperature [33| ; when an impurity scatters from a set 
of different fixed impurities with scattering resulting in a 
Fermi distribution in momentum space [34j | ; and where 
an impurity scatters on phonons |35j . 

Turning this on its head, the precise way in which the 
current of the impurity deviates from the Esaki- Tsu de- 
pendence reveals information about the environment and 
impurity-environment interaction. For example, the abil- 
ity to excite optical phonons in the environment at a 
single frequency result in resonances in the current of an 
impurity at certain tilts dependent on this frequency [35| . 
The interaction of impurities with a supcrfluid in a lat- 
tice, which has a divergence in the density of phonon 
states at the band edge, also causes peaks at forcings 
related to fractions of the bandwidth j8[. 

We now consider in which way properties of the en- 
vironment manifest themselves in the current of the im- 
purities, approaching the problem using both the semi- 
classical and Wannier- Stark pictures. 



Scattering in momentum space 




tJJh 



100 



(b) 0.5 




rA/h 



FIG. 2: (Color online) The motion of a single impurity in a 
tilted lattice with and without a scattering process, (a) The im- 
purity decoupled from any environment will coherently Bloch 
oscillate at frequency A/h but its expected position will not 
drift [1, H3|. Here we plot the evolution in time t for just over 
one and a half Bloch oscillations. The impurity is initially 
localized in a Wannier state at the middle site of the 101-site 
lattice. We set the tilt to be A/J a = 0.1 and (n a ,i) is the 
impurity density at lattice site i. (b) Introducing a scattering 
mechanism under the relaxation-time approximation results 
in the Esaki- Tsu dependence of drift on tilt in Eq. Q. We 
plot this dependence here. 



obeys the Boltzmann equation [2 
dg(k) 



dk 



ha 

a" 



dq[W(q,k)g(q)-W(k,q)g(k)], (4) 



which is valid for homogeneous systems of impurities in 
the steady state and where W(q, k) is the rate of inco- 
herent scattering between quasi- momenta q and k. The 
average velocity of each impurity is then the expected 
value of the group velocity of wave-packets distributed 
according to g(k), 



dkg(k)v(k), 



(5) 



First, we use the semi-classical approach where the dis- 
tribution of impurities in quasi- momentum space g{k) 



where the group velocity for a particular quasi- 
momcntum is defined in Eq. @. For small enough 



4 



tilts the quasi-momentum will be a very slowly changing 
quantity and hence for weak interactions with the envi- 
ronment we may approximate the incoherent scattering 
rates using Fermi's golden rule as 

2tt 



W(q,k) =— \(k\(n b \H ab \q)\O b ) 



x 5 (E nb — E 0b — 2J a [cosqa — cosfca]) . 

Here we have assumed the environment to be sufficiently 
large and cold so that before scattering we may take it 
to be in the ground state | 0&) of energy Eo b . The result 
of the interaction of the impurity with the environment, 
for which the Hamiltonian is given by H abl is to scatter 
the environment to some excited state | n b ) with energy 

En b - 

Particularly important for our system is the case of 
impurities coupled to a many-body environment, here 
the Bosc gas, via a density-density interaction. For low 
energies, the interaction takes the explicit form 



H a b = U ab ^ a\a,ib\bi 



(G) 



site i is normalized to Yli Pi(t) 
average velocity then reads 



1. The corresponding 



<«>=E' 



■ dP«(t) 
1 dt ■ 



Note that in an infinite homogeneous system the average 
velocity would be (v) = J2j a ti ~ i)Wij. Again, for weak 
coupling to the environment the transition rates Wij be- 
tween different Wannicr-Stark states are given by Fermi's 
golden rule as 



^■ = T^I (jK "' b| ^ ab|l)l ° b) 

n b 

x5(E nb -E 0b + [i-j]A). 



(8) 



In this picture negative differential current arises due to 
the increasing localization of the Wannier-Stark states: 
as A increases, the width of the Wannier-Stark states 
decreases as 1/A resulting in a suppression of their over- 
laps and thus the motion of the impurity. 

Inserting the interaction Hamiltonian in Eq. ^ into 
Eq. ([5]) we arrive at the environment correlation function 



where U ab is its characteristic energy and b\ (hi) is the 
creation (annihilation) operator for a particle in the en- 
vironment at lattice site i. The corresponding quasi- 
momentum scattering rates are given by [361 ] 



2nU' 



W(q,k) = ^^S(q 



fe, — 2 J a [cosqa — cosfca]), (7) 



where L is the length of the system, S(q, lu) is the zero- 
temperature dynamic structure factor of the environment 

S(q,Lo) = J2\(n b \p q \0 b )\ 2 6(E nb -E Qb -fkj), 



and p q is the Fourier transform of the particle density of 
the environment. Hence the impurity current and, more 
directly, the momentum distribution at small tilts are 
governed by the dynamic structure factor of the environ- 
ment. 



2. Transitions between Wannier-Stark states 

At larger lattice tilts it is convenient to work in the 
basis of Wannier-Stark states and to consider the inco- 
herent transitions between these. The dynamics of the 
impurities is described by the master equation 



dPj(t) 
dt 



with the incoherent hopping rates W^. The probability 
Pi that the impurity occupies Wannier-Stark states at 



EME Ji-jWt-iWbtbel b ) 



xS(E nb -E 0b + [i-j]A). 



. (9) 

which determines the current at larger tilts. From this we 
immediately see two effects that we expect to be visible 
in the impurity current. First, consider the case where 
the environment has a gap G in its excitation spectrum 
such that E nb — Eo b > G. This means that the delta 
function in Eq. © will be zero unless A > G/(j — i) and 
hence the contributions to the current due to hopping 
between Wannier-Stark states shifted by I sites will only 
be made at tilts above A = G/l. Second, if the excita- 
tion spectrum of the environment has a finite bandwidth 
such that there are no states above an energy B, i.e., 
E nb — Eo b < B, then we expect ^-site contributions to 
the current to be suppressed at tilts above A = Bjl. 
Hence the deviation of the dissipation process from the 
relaxation-time approximation results in the dependence 
of the impurity current on A exhibiting sharp features 
that clearly reveal information about the environment; 
we can infer these environmental properties by observing 
the impurities. 



C. Properties of the Bose gas environment 

Our specific environment consists of a Bose gas trapped 
in a horizontal (non-tilted) optical lattice in which the 
bosons are confined to move in the same single dimen- 
sion as the impurities. As for the impurities we consider 
temperatures to be low enough such that only the low- 
est Bloch band need be considered. In such a case the 
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bosons are described by the Bose-Hubbard model |3j 



where the operator b\ (5,) creates (annihilates) a boson 
in a Wannier state localized at site i and the lattice pa- 
rameter a is the same as that for the impurities. Jb and 
Ub are parameters that determine the hopping between 
neighboring sites and the on-site interaction, respectively, 
and can be tuned experimentally by adjusting the laser 
parameters or using Feshbach resonances [H, [3?} ■ 

Let us briefly review some properties of the Bose gas 
which we will use to interpret our findings. In the super- 
fluid regime (Ub/ Jb <C 1) the Bose gas supports excita- 
tions in the form of Bogoliubov phonons, whose excita- 
tion spectrum is given by 



huj k = v /e fc (e fc + 2U b n b ), 
Ek = 2Jb(l — coska), 



(10a) 
(10b) 



where rib is the ratio of bosons to lattice sites. From 
this we find that the density of states (hdujk/dk)^ 1 di- 
verges at the upper edge of the band with bandwidth 
B = AJb\Jl + Ubrib/2Jb. In the opposite limit Ub/Jb — > 
co the bosons map to the same number of non- interacting 
identical fermions with energy spectrum e^. Hence in this 
regime the Bose gas supports particlc-hole excitations for 
nb < 1 within the single band of bandwidth B = AJb- 

For the commensurately filled case of n& = 1 an in- 
creasing interaction strength Ub/Jb takes the gas through 
a continuous superfluid to Mott insulator transition, 
which occurs at Ub/Jb = 3.37 [HI, S3 m the thermo- 
dynamic limit. A principal signature of this transition 
is the appearance of a gapped excitation spectrum sup- 
pressing the low-energy response. Deep in the Mott in- 
sulator regime (Ub/Jb ^ 1) the gap in the excitation 
spectrum is G ~ Ub- 



D. Detailed model and measuring procedure 

To summarize, the model for our system is given by 
the impurity Hamiltonian H a , the Bose- Hubbard model 
Hb and the density-density interaction H a b, as shown 
schematically in Fig. [T] More precisely, we consider an 
Af-site system with box boundary conditions or a suit- 
ably flat bottomed potential (see |4l[ for an example of a 
box trap realized in an experiment). As in Q we expect 
our system to share its bulk behavior with the experi- 
mentally important case of an additional harmonic trap, 
provided that the trap is sufficiently shallow. A realiza- 
tion of the model may also draw on recent experimental 
successes in trapping atoms in species-specific optical lat- 
tices mi sup. 

The system is identical to that considered in @ which 
focussed on a superfluid environment (Ub/Jb <C 1), 
whereby the impurity and Bose gas together mimicked 



the clcctron-phonon interaction. However, the control 
and flexibility of an optical lattice setup allows for the 
creation of a strongly- interacting Bose gas (Ub/Jb 3> 1), 
the regime we explore in this article. 

To measure the impurity current we propose the follow- 
ing procedure: The N bosons are cooled to their ground 
state with the impurities fixed in highly localized states 
at their initial locations (for example, using a tight opti- 
cal trap). The system reaches equilibrium and the pres- 
ence of the impurities inevitably leads to density varia- 
tions in the initial state of the bosons as compared to the 
ground state of the Bose-Hubbard Hamiltonian. At time 
t = the impurities arc released, their lattice is tilted and 
the system is left to evolve. Later, a time of flight mea- 
surement or in situ image 0, EH is taken [52[ • Note 
that alternative preparations of the initial state we con- 
sidered, e.g., preparing the bosons in equilibrium without 
the impurities and then adding the impurities at t = 0, 
produced the same qualitative results despite this being 
a more energetic initial state. 

For our numerical simulations we consider the com- 
plete two-species system governed by the Hamiltonian 
H a + Hb + H a b- The simulations of the full many-body 
quantum dynamics at zero temperature are performed 
using the TEBD algorithm. This is a matrix prod- 
uct based method extending the powerful density matrix 
renormalization group to real time evolution. We present 
an overview of the method in Appendix [Al along with the 
algorithm parameters used. 

We consider a single impurity initially localized in a 
Wannier state at the center of the lattice. We evolve the 
system up to t s i m ss Mh/AJb, the time it approximately 
takes the disturbances in the bath caused by the impurity 
to reach the boundaries of the system, thus minimizing 
finite size effects. This time limit is determined by the 
maximum group velocity v g of the bosons in the band 
which can be derived in both the superfluid Ub/Jb 
and strongly- interacting Ub/Jb — > oo limits by differenti- 
ating hwk and Ek [Eqns. (jTDJ)] , respectively, with respect 
to k. The time taken for disturbances in the Bose gas 
to travel a distance Ma/2 to the edge of the system is 
then Ma/2v g from which the previous expression for t s i m 
follows. We choose to hold the interaction between the 
impurity and the bosons constant at U a b/Jb = 0.5 and 
the impurity hopping at J a /Jb = 0.6. These values are 
chosen to reflect the behavior we observed over a wide 
range. Notably they are outside the perturbative regime 
used in Sec. Ql] yet still lead to the predicted effects. 



III. SIMULATIONS OF IMPURITY 
PROPAGATION 

As a demonstration of the system and numerical 
method we show the evolution of both the impurity and 
bosons in Fig. [3] for an intermediate boson-boson interac- 
tion strength Ub/Jb — 2. The impurity undergoes Bloch 
oscillations as it would if the bosons were not present — 
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FIG. 3: (Color online) The transport of an impurity through 
a bosonic bath. An example of the evolution of the impurity 
and boson densities for the case A/Jj = 0.3, Ub/Jb = 2, 
M = 101 and N — 50. The color map in the xy-pl&ne shows 
the expected impurity density at each lattice site (h a ,i), and 
similarly for the surface plot and boson density 

see Fig. [2ja) — however, the non-zero impurity-boson 
interactions cause the impurity to dissipate energy into 
the Bose gas and drift down the lattice. The same inter- 
actions also create density fluctuations in the Bose gas, 
which spread out and reach the system boundary near 
the end of the simulation. 

Examining the system for a range of tilts A allows us 
to find the dependence of impurity displacement (x a ) = 
a{io — i)(fi>b,i) on A, where zo is the site at which the 
impurity was initially localized. Using this we can calcu- 
late the average impurity velocity (v) = {x a )/t s i m . From 
detailed examination of our simulations we find that the 
displacement does not increase just linearly in time but 
oscillates slightly at the frequency of the Bloch oscilla- 
tions (as would be observed in an experiment). For all 
but the lowest tilts the durations of evolutions considered 
is sufficient that by i s i m the displacement has averaged 
over enough Bloch oscillations that (v) gives a good rep- 
resentation of the long-time average drift. 

A. Superfluid and Tonks-Girardeau regimes 

Our starting point is the propagation of the impurity 
through a Bose gas in the superfluid regime with incom- 
mensurate filling. A typical dependence of the displace- 
ment of the impurity (x a ) on the lattice tilt A is given 
by the red circles in Fig. HJa). This dependence has the 
characteristic features of the Esaki-Tsu result, shown in 
Fig. [5fb), namely an ohmic region followed by negative 
differential conductance. 

However, a phenomenon that is specific to the case 
where the bosons feel a lattice is the occurrence of res- 
onances in the vicinity of the lattice tilts A = AJi,/i for 



(a) A/J b (for MI) 




FIG. 4: (Color online) Determining the bandwidth and gap 
from the dependence of displacement on tilt, (a) The depen- 
dence of the displacement of the impurity (x a ) at i s im as a 
function of tilt A for three regimes of the Bose gas: the su- 
perfluid regime U b /J b = 0.1 with N = 50, M = 101 (SF 
- displacement scaled down by a factor of 0.1); the Tonks- 
Girardeau limit Ub/Jb -> oo with N = 150, M = 151 (TG); 
and the Mott insulating regime U b /J b = 10 with N — M — 
101 (MI - displacement scaled up by a factor of 100). (b) 
The same quantity is plotted for several bosonic systems in 
the Tonks-Girardeau limit, demonstrating that the qualita- 
tive shape and the positions of the bandwidth suppressions 
are independent of the filling (so long as it is incommensu- 
rate) of the Bose gas. 



integer £ > 1. As shown in [8( we expect resonances at 
tilts A = (4Jb/£)y/l + f/&nb/2Jb corresponding to a di- 
vergence in the density of phonon states at the upper 
band edge. Moreover, the suppression of the propaga- 
tion of the impurity for A > AJ^/t reflects the sudden 
inability of the impurity to move £ sites down the lattice 
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by dissipating the energy £A into the Bose gas through a 
single phonon process. Hence the superfluid Bose gas is 
a prime example of how the value of the bandwidth of its 
excitation spectrum B manifests itself in the dependence 
of the current of the impurity. Note that the small fluc- 
tuations in the dependence of drift on tilt at small tilts in 
Fig. Hta) are a result of the finite time of our simulations 
and would be smoothed out by considering longer times. 
The same would be true in an experimental realization. 

One expects to see the bandwidth dependent features 
even with strong boson-boson interactions; the reasoning 
used to predict current drop offs after A = B/i in SecHTl 
assumes the bosonic system is susceptible to low energy 
excitations, which is the case for incommensurately filled 
systems of any interaction strength [38j. Indeed, we have 
simulated M = 25 site systems over a range of boson- 
boson interaction strengths Ub/Jb (not shown here) and 
found no qualitative change in the behavior of the mo- 
tion. However, as shown in Fig. [5j the average displace- 
ment of the impurity decreases with increasing interac- 
tion strength. 

In the extreme case of infinitely strong interactions 
Ub/Jb —> oo, i.e., in the Tonks-Girardeau limit, the situa- 
tion simplifies considerably. The bosons map to identical 
non- interacting fermions with bandwidth B = AJb, and 
so the creation of particle-holes at the Fermi surface pro- 
vides low-energy excitations. Accordingly, we expect a 
large impurity displacement at low lattice tilts and sharp 
drops in propagation for tilts above integer divisions of 
the bandwidth 4Jb/£, in agreement with Fig. B|a). 

We find the same behavior in the Tonks-Girardeau 
limit for many different bosonic filling factors, as shown 
in Fig. 0Jb). The half- filled bosonic system gives rise to 
a peak impurity current that is an order of magnitude 
above that for the cases where the bosonic lattice is filled 
by a single or M — 1 bosons. This can be attributed to 
there being a greater number of excitations available in 
the half-filled bosonic system; using the fcrmionic map- 
ping we find there arc of approximately M 2 /4 single 
quasi-particle excitations available in the half-filled case 
as compared to M — 1 when there is only one boson or 
M — 1 present. This reasoning suggests that the peak 
displacement may be used as a measure of the filling of 
the bosonic system. 



B. Superfluid to Mott insulator transition 

We now turn to the case of a commensurately filled 
Bose gas and the continuous superfluid to Mott insula- 
tor transition. Since motion of the impurity though the 
superfluid at low lattice tilts is a result of the ability of 
the bosons to accept low-energy excitations we expect 
propagation to be highly suppressed in the Mott insula- 
tor phase which has a gap in the excitation spectrum of 
G~U b . 

The effect of this signature on the propagation can be 
seen in Fig. HJa) , where the displacement has a peak at 
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FIG. 5: (Color online) The suppression of average drift in the 
strongly-interacting regimes. Here we show the peak displace- 
ment (maximized over all tilts) as a function of interaction 
strength Ub/Jb for commensurate and incommensurate fill- 
ings of a 101-site bosonic lattice. 



A « Ub- There is also a second peak at A ss Ub/2 result- 
ing from second-order processes involving two lattice-site 
jumps, similar to what is seen using Bragg spectroscopy 
in [19| and tilting the bosonic lattice in [21j. As in those 
cases the non- zero Jb results in a broadening of the peaks. 
This behavior confirms the predictions in Sec. |H]for the 
transport of an impurity in a gapped environment. This 
drastic change from the behavior predicted by Esaki and 
Tsu also tells us that the relaxation-time approximation 
is not at all valid for the interaction of an impurity with 
a Mott insulator. 

For small tilts the motion of the impurity caused by 
low-frequency excitations is completely suppressed; at 
large enough Ub/Jb the ground state is well approximated 
by a unit-filled Fock state. In such a regime the impurity 
Bloch oscillates as it has no means of dissipating energy 
into the Bose gas, hence there is no net current across 
the lattice. As a result we get a startlingly clear signa- 
ture of the superfluid to Mott insulator transition in our 
dynamics in the form of a quench of more than 4 orders 
of magnitude in the propagation of the impurity. This is 
shown in Fig. [5] 

As well as the ability to probe the superfluid to Mott 
insulator transition using the impurity, this quench pro- 
vides a non-destructive way of establishing whether the 
bosonic system is commensurately filled or not. If our 
system is strongly interacting Ub/Jb ^ 4 there are several 
orders of magnitude separating the currents for an im- 
purity in the commensurate and incommensurate cases, 
shown in Fig. [5] Many implementations of quantum in- 
formation processing involve creating a commensurately 
filled system. Using an impurity as a probe is one way 
in which this commensurability could be checked non- 
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destructively. 




FIG. 6: (Color online) Quasi-momentum distribution of the 
impurity after traveling through the Bose gas. Here A/Jb = 
0.063, h a . p is the number operator for an impurity quasi- 
momentum of 2np/aM and M = 101. (a) The momentum 
distribution of the impurity at t s i m for a range of interaction 
strengths Ub / Jb and incommensurate filling TV = 50. (b) The 
same as (a) but for commensurate filling TV = 101. 



C. Quasi-momentum distribution of the impurity 

Alternatively, one can gain information about the im- 
purity by looking at its quasi-momentum distribution, 
which is directly accessible in time-of-flight experiments. 
With the impurity prepared in a Wannier state each 
quasi-momentum state is equally occupied. However, 
we expect interactions with the Bose gas to alter this 
momentum distribution, resulting in the drift seen in 
Sec. IIIII The current (v) can be calculated from the dis- 
tribution in the semi-classical picture using Eq. (|5|) . 

As can be seen in Figs. Uta) and [Bfb) the result of 
interactions with the Bose gas is the enhancement of the 
distribution near the centre of the Brillouin zone. More- 
over, the interaction gives rise to an asymmetry with the 
distribution skewed towards negative quasi-momentum, 
necessary for the impurity to have a current. 

Our simulations show that a well-defined quasi- 
momentum emerges in the superfluid regime by the time 



t s ; m , supporting the accuracy of the semi-classical ap- 
proximation. However, for stronger bosonic interaction 
strengths the scattering results in a wide distribution 
with a negative quasi-momentum shoulder. By Ub/Jb 
of the order 10 the distribution already approximates its 
Ub/Jb — > oo value. In the case of the Mott insulating 
regime in Fig. [B^b) it is unchanged from its initially flat 
distribution. Once again, the momentum distribution of 
the impurity provides a signature of the transition (cross- 
over) between the superfluid and Mott insulator (Tonks- 
Girardeau) regimes. 

The difference in steady-state quasi-momentum distri- 
bution is related to the structure factor of the Bose gas, 
c.f. Eqns. (HJ) and ([7]), but a quantitive comparison with 
our simulations is difficult since in general they do not 
reach the steady-state. Since this limitation is due to the 
finite time of our simulations it would also restrict an 
experimental realization. 



D. Reduced density matrix of the impurity 

So far we have focused on the evolution of the spatial 
and quasi-momentum distributions of the impurity, how- 
ever more information is contained in its reduced density 
matrix p a . First, evaluating the elements of the den- 
sity matrix in the Wannier basis may tell us more 
about the effect of boson-impurity interactions on the 
impurity. Another quantity that can be calculated from 
the reduced density matrix is the purity Tr(/3^), which 
takes the maximal value 1 if p a describes a pure state 
and decreases for increasingly mixed states. The initial 
state of the total system is a product of pure states for 
the impurity and bosons, and since we consider coherent 
evolution it remains pure. However, through interactions 
of the impurity with the Bose gas the two subsystems 
are entangled resulting in a loss of purity of the impu- 
rity. Our bosonic environment is sufficiently large that 
for most parameter values revivals are not observed and 
so we use purity to quantify the decoherence of the im- 
purity. 

To demonstrate the behavior of the density matrix p a 
we have plotted in Figs. EJa) and[7{b) the absolute values 
of its Wannier basis elements at t s ; m for both intermedi- 
ate and weak bosonic interaction strengths, respectively. 
In Fig.[7Ja) we observe a symmetric coherent interference 
pattern in the off-diagonal elements as the impurity un- 
dergoes Bloch oscillations. Additionally, density-density 
interactions with the bosons result in a slight decay of 
the off-diagonal coherences and an asymmetric distribu- 
tion along the diagonal with a higher occupancy of the 
lower energy lattice sites. This decoherence is much more 
pronounced for a superfluid environment, as shown in 
Fig. Efb), and the resulting density matrix is nearly di- 
agonal. 

To make this relationship between the decoherence of 
the impurity and the internal interaction strength of the 
environment clearer, Fig. [7f c ) shows the purity of the 
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FIG. 7: (Color online) Decoherence of an impurity traveling 
through an incommensurately filled system N = 50, M = 101. 
(a) Absolute values of the elements of the reduced density 
matrix for the impurity |(oJa,-)| with Ub/Jb = 2 and A/Jf, = 
0.1, calculated at t s i m . (b) The same as (a) but for Ub/Jb — 
0.1. (c) Also at t s i m we have plotted the purity of the density 
matrix Tr(p^) as a function of tilt and interaction strength. 



reduced density matrix of the impurity as a function of 
both tilt A and interaction strength Ub/Jb- As can be 
seen in the figure, the coherence of the impurity is a clear 
signature of the cross-over between the superfiuid and 
Tonks-Girardcau regimes. Moreover, the values of Ub/Jb 
at which the most significant changes occur is around the 



same values where the transition occurs in a commen- 
surate Bose gas. We also calculated the entanglement 
entropy between the impurity and Bose gas (not shown 
here) which showed very similar behavior to purity. If 
the Bose gas is commensurately filled the transition be- 
tween weak and strong interaction strengths makes an 
even clearer impact on the purity as deep in the Mott 
insulator phase the impurity does not become at all en- 
tangled with the Bose gas. 

Finally, comparing Figs. [7|c) and [3jd) we find that 
the regimes of high transport coincide with those of the 
smallest purity; the very same scattering processes that 
lead to the propagation of the impurity along the lattice 
also cause its decoherence. Therefore the propagation 
of the impurity provides a controllable example of noise- 
enhanced transport. Noise-enhanced transport is a topic 
of interest for those searching for quantum effects at room 
temperature which explain the remarkable efficiency of 
some biological systems [46[ . 

IV. CONCLUSIONS 

We performed fully quantum many-body simulations 
of large (100 or more sites) two-species systems using the 
TEBD algorithm, simulating the one-dimensional motion 
of an impurity confined to a tilted optical lattice through 
a Bose lattice gas. We discussed the Esaki-Tsu depen- 
dence of the current of the impurity on tilt and provided 
evidence that this shape is qualitatively obeyed by an 
impurity traveling through an incommensurately filled 
bosonic system of any interaction strength including the 
Tonks-Girardcau limit. Contrasting this a very different 
dependence was found in the Mott insulating regime. 

We discussed how deviations to the Esaki-Tsu shape 
depend on properties of the Bose gas. Information about 
the gas can then be extracted without needing to mea- 
sure the bosons themselves, i.e., non-destructively. In 
particular we discussed how the impurity could be used 
to probe the commcnsuratcncss of the system, and obtain 
the gap and bandwidth of its excitation spectrum. 

Following this we analyzed the transport of the im- 
purity by considering its momentum distribution, show- 
ing that the semi-classical picture of an impurity of well- 
defined quasi-momentum emerges after a short time for 
bosons in the superfiuid regime. We have also found a re- 
lationship between the decoherence and transport of the 
impurity; our setup provides an experimentally accessible 
prototype of noise enhanced transport. 

Appendix A: Time-evolving block decimation 

The TEBD algorithm efficiently simulates both the 
unitary and imaginary time-evolution (according to 
c -iHt anc j Q-Ht rcS p C ctivcly) of quantum systems com- 
prising a one-dimensional lattice of subsystems, each with 
a finite number of configurations, where H is composed 
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of nearest neighbor terms [lj, [l5[ . It does this by storing 
the state of a system as a matrix product state (MPS) of 
dimension \. In general, for an exact representation of 
the state vector the dimension \ and hence the resources 
required for the simulation must grow exponentially with 
the size of the lattice. However for ground and low-lying 
excited states of one-dimensional systems near-exact ac- 
curacy can be obtained using a much smaller v that does 
not grow exponentially with system size 4?], 48 1. With 
this being the case and using a Suzuki- Trotter decompo- 
sition of the operator e~ lHSt (without the i for imaginary 
time-evolution) [4j| the evolution may be simulated by 
the application of a polynomial number of two-site gates, 
each followed by a re-compression into MPS form, which 
can be done efficiently. For details, we refer the reader 

to m M, M M- 

Our system, described in Sec. lU is of this type pro- 
vided we restrict the occupancy of the bosonic lattice 
to a finite number. To realize our scheme we first used 
TEBD to find the ground state of the bosonic part of 
the system, with an on-site potential U a b added at each 
site corresponding to an impurity. The initial state of 
the whole system was then the tensor product of this 
bosonic state with the initial state of the impurity, which 
consists of Fock states with occupancies or 1. Next, 



TEBD was used to simulate the evolution of this state 
under the full Hamiltonian H = H a + Hb-\-H a b. Expecta- 
tion values such as bosonic and impurity densities were 
easily extracted from the numerics due to the efficient 
contractibility of an MPS. 

For our simulations we used a maximum boson lattice- 
site occupancy of 4 and the largest MPS dimension we 
used was \ = 120. Our timestep was St = 5 x 1Q~ 3 H/Jb- 
These parameters were sufficient for the time-evolution 
considered here; we found that increasing \ an d decreas- 
ing St resulted in no significant changes to the observables 
computed. 
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